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THE OBERBECK-BOUSSINESQ APPROXIMATION IN CRITICAL 

SPACES 

RAPHAEL DANCHIN AND LINGBING HE 

Abstract. In this paper we study the vahdity of the so-called Oberbeck-Boussinesq ap- 
proximation for compressible viscous perfect gases in the whole three-dimensional space. 
Both the cases of fluids with positive heat conductivity and zero conductivity are consid- 
ered. For small perturbations of a constant equilibrium, we establish the global existence of 
(~^ • unique strong solutions in a critical regularity functional framework. Next, taking advantage 

^S] I of Strichartz estimates for the associated system of acoustic waves, and of uniform estimates 

with respect to the Mach number, we obtain all-time convergence to the Boussinesq system 
O ,' with a explicit decay rate. 

1. Introduction 

^ . This work aims at giving a mathematical justification of the Oberbeck-Boussinesq approx- 

imation that is commonly used to model stratified fluids such as e.g. atmosphere or oceans. 
One of the characteristics of the this approximation is that, although the primitive system 
2 ■ is the full compressible Navier-Stokes system, the limit equations are incompressible, and 

the density is a constant. In fact, the velocity field just convects an active scalar creating 
buoyancy force, proportional to the discrepancy between the temperature and its equilibrium. 

, 1^ I 1.1. Formal derivation. The starting point of our analysis is the full Navier-Stokes system 

r-^ ■ for compressible viscous fluids, namely 

( dtp + div (pu) = 0, 



^ 



o 

en 



1 1 

dt{pu) -\- div {pu (g) u) - divr H ^VP = -^ pW, 

Ma Pr 

^ dt{ps) -\- div (ps) -\- div {q/T) = a. 



Above p = p{t,x) G R^, u = u{t,x) G M and T = T{t,x) G M"*" stand for the density, 

S^ I velocity field and temperature, respectively. The scalar function V stands for some (given) 

H ■ external potential (e.g. the gravity potential). We concentrate on the study of the evolution 

- - - toward the future in the whole space M (hence the time variable t belongs to M"*" and the 

space variable x, to M ). 

In the Newtonian case that we shall consider, the stress tensor r is given by 

r = piVu + Du) -\- Adiv u Id. 

For simplicity, the viscosity coefficients A and p are assumed to be constant. As we only 
consider viscous fluids, those two coefficients satisfy 

/_i > and V := \ + 2p> Q. 

This ensures ellipticity for the second order operator A := pA -I- (A -|- ^)Vdiv . 
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The heat flux q is equal to —hc^T for some constant conductivity coefficient k > 0. The 
pressure P, the internal energy e and the specific entropy s are related to p and T through 
the Gibbs relation 

Tds = de + Pd{l/p). 

We focus on perfect gases, namely we assume that for some a > and 6 > 0, P = apT 
and e = bT. After rescaling, it is non restrictive to take a = 6 = 1. 

Finally, in the velocity equation, the Mach number Ma and the Froude number Fr are 
two dimensionless small parameters accounting for the compressibility and the stratification 
of the fluid. Formally, Oberbeck-Boussinesq approximation is obtained in the asymptotics 
e^Oif 

Ma = e and Fr = \/e, 

an assumption that we shall make from now on. 

Gathering all the above assumptions over the coefficients and state laws, we end up with 
the following system (with exponents e emphasizing the dependency with respect to e): 

' dtp' + div {p'u') = 0, 



(1.1) 



dtip'u") + div {p^u" (g) u" 
dtip'T') + div {u' p'T') 



VP^ 1 

- pAu^ - (A + pWdivu'' H ^ = -p^Vy^ 

e^ e 

kAT' = e^[2p\Du^\'^ + A(divn^)2]. 



Let us first provide a formal derivation of the Oberbeck-Boussinesq approximation in the 
case where the heat conductivity k is positive. We want to consider so-called ill-prepared 
data of the form Pq = 1 + ea^, Uq and Tq = 1 + eOq where (oq, ^^Q' ^o) ^^^ bounded in a sense 
that will be specified later on. Setting p'^ = 1 + ea^ and T^ = 1 + eO'^, we get the following 
governing equations for (a^,ti'^,0^): 

„ , div u'^ 
dta"" + 



(1.2) 



dtu'' + u" ■ Vu" 



-div (a^u^ 
1 + ea^ 



+ 



V(a= + r + ea^ 



1 



£(1 + ea^ 



W 



dte' + 



divu*" 



-I- div I 



3'u') 



kA9^ 



[2p\Du''\'^ + Xidivii: 



e\2i 



e ■ ' 1 + ea*^ 1 -I- ea'^ 

In order to handle the singular potential term in the r.h.s. of the velocity equation, it is 
usual to work with the modified deviation of density b^ := a^ — V^ . We get 

divn*" 



(1.3) 



dtb' + u' ■ Vb' + 



dtu^ + u^ ■ Vu^ - Au^ + 
divti^ 



-dtV - div {V'u") - 6Mivn^ 



dtO' 



vr + 



v(6^ + r 

e 

- KAr = 



l + ea^ 



Va^ 



ea' 



l + ea" 



Au^, 



e\2] 



1 + ea^ 



[2/i|Z)u^|^ + A(divn^)^] 



ea 



l + ea^ 



Ar-rdivu^ 



which may formally written as follows: 

b^ \ . / div 

+ - I V 

div 





0(1). 
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The notation 0(1) designates terms that are expected to be bounded uniformly with respect 
to e. 



As a consequence of our considering ill-prepared data, the first order time derivatives are 
likely to blow-up like 1/e for e going to 0. At the 'physical' level, this means that highly 
oscillating acoustic waves may propagate in the fluid. 

In order to better understand the action of those singular terms, we may first look at 
the kernel KerL of the 5x5 first order antisymmetric differential matrix operator L above. 
The basic idea is that modes that are in Ker L will not be affected, while modes that are in 
(KerL) may experience wild oscillations. A straightforward computation shows that 



KerL 

(KerL)- 



|(6,u,6l) :divn = and V(6 + 6l) = o|, 
|(6,u,6l) :curln = and V(6-6l) = o|- 



Hence it is natural to look more closely at the equations satisfied by {q^ , Qu^) and {Q'^jVu^) 
where V and Q stand for the orthogonal projectors over divergence-free and curl- free vector 
fields, respectively, and 

e^ + b"" „. 9^-b^ 



V2 ' 



e^ 



^/2 



As L is antisymmetric, we expect the oscillating components of the solution, namely Qu^ and 
q'^ to be dispersed whereas L will have no effect on Vu^ and 0"^. Let us be more accurate: we 
see that {q^, Qu'^) satisfies 



(1.4) 



dtq' + — div Qu' = -div (q'u') - — ( dtV + div (V'u') + k ^^' 
£ 2 V 1 + ea^ 



+ 



V2 



dtQu' + ^Vq' = Q 



l + ea^ 



Va" 



2 l + ea^ 



[2/i|L»'u^|^ + A(divu' 



e\2i 



1-Fea^ 



u' ■ Vu' 



whereas (G^,Pn^) fulfills 



(1.5) 



ae^ + pM^-vG^--Ae 



/2 



Ag" 



V2, 



K sa 



-Ar + 



V2 



e\2^ 



2^|Z)n^|^ + A(divu^) 



2 2 l + ea^ 2 1 + ea^ ' 

dtVu^ - fiAVu'' + V{Vu^ ■ VVu^) + P(rVa") = -V{u^ ■ VQu^ + Qu^ ■ VVu^) 



V 



ea^ 



l + ea'^ 



Aun+v 



ea^ 



a^) 



l + ea'' 



Va' 



If we assume the solution {b^ ,u^, 6^) and the data to be bounded independently of e then 
the right-hand side of (|1.4p is bounded, too. Hence, owing to the antisymmetric (and nonde- 
generate) structure of the left-hand side of (jl.4p . one may expect {q^, Qu^) to tend weakly to 
0. We shall see later on that in the whole space setting that is here considered, it is possible 
to get strong convergence (for suitable negative Besov norms), with an explicit rate. 
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In order to find out what the hmit system for (|1.5p is, let us observe that 

\/2p(rva^) = r{e'vv')+r{q'W) + V2V{e'vw) 

= P(G^Vy^) + V{q'VV') + \/2P((r + ¥)Vh^) 

(1.6) = V{Q'VV^) + V{q'VV^) + 2V{q'Vb^). 

Because q^ tends to 0, we expect that 

V2V{e^Va^) -V{e'VV') ^ for e going to 0. 

Hence, if we assume in addition that V^ — )• V, Vuq — ;■ vq and 0q — )■ Gq, then (0'^,'Pu^) should 
tend to the solution (B,f) to the following Boussinesq system: 

^/2, 



(1.7) 



dtQ + v-VQ- -Ae 



-(9i + r;-V)y, 
V2, 



ati; + 1; • Vt) - ^Au + vn = --^evy, 



divf = 0, 



JG,t;)|t=o = (eo,^o). 

Setting = — v2/2 y, and changing VII accordingly, we see that this system is equivalent 
to the following one, which is commonly used: 



(1.^ 



dtQ + v-VQ- -^AO = ^kAF, 

dtv + vVv- ^Av + vn = --^ovy. 



divf = 0. 



Note that although the density is constant in the limit system, it comes into play in the 
buoyancy force where it is related to the temperature and the potential. 

We end this paragraph with a formal derivation in the case k = 0. It turns out to be easier 
to work with the pressure rather than with the temperature. We thus set p^ = 1 + ea^ and 
P^ = p^T^ = 1 + e{n^ + V), and obtain that 

divti*^ 



(1.9) 



dta' + 



dtv!" + n^ • Vu" 



-div(a'^ti^), 
+ 



VW 



1 + ea^ e(l + ea=) 1 + ea< 



■W 



dtn' + 



divn^ 



£\2l 



+ div {n^u") = e[2p\Du''\^ + A(divu^)^] - dtV^ - div {V'u'). 



Setting 0^ := a'^ — IZ"^ — V^ , we thus get 

(dt@' + div (O^u") = -e[2p\Du^\^ + A(divu=)2], 



(1.10) <^ 



dtVu^+V{u^ ■ Vu') - pAVu' 



dtQu^ + Qiu^-Vu") - z^AQn^-f- 



-V 



VTZ' 



ea^ 



l + ea= 



Aun+v 



l + ea= 



-Q 



ea' 



l+ea^ 



Au^+Q 



ViV'+W 



v{v^+n^ 



l+ea^ 



dtn' + 



divn^ 



£\2^ 



+ div(7^"n^) = e[2p\Du'\^ + A(divM^)^] - dtV - diw{V'u'). 



As before, owing to the first order antisymmetric terms, we expect {Qu^ ,1Z^) to go to 0. 
Concerning (0^,^X1^), we notice that 
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Therefore the hmit system for (B^,'Pn'^) reads 

' 9tG + u • VG = 0, 
(1.11) J dtv + v-Vv- ^Av + Vn = QVV, 

divw = 0. 

Note that in contrast with (jl.7p . this system is not fuhy parabohc. 

1.2. Some related ■works. There is an important hterature dedicated to the Hmit system, 
that is the Oberbeck-Boussinesq equations (II. Th . (II. Sh and (II. lip , under various hypotheses 
over the coefficients k and //, and the potential V (although the most common assumption is 
that V = X3). Loosely speaking the classical results concerning the existence issue are (see 
e.g. [8l[9l[16] and the references therein): 

• Dimension 2: Global existence of strong solutions if (//,«) 7^ (0,0). 

• Dimension 3 with fj, ^ 0: Global weak solutions and local strong solutions (which 
become global if the data are small). 

• Dimension 3 with fi = : only local-in-time strong solutions are available. 

In contrast, although the Oberbeck-Boussinesq approximation is commonly used in geo- 
physics (see e.g. the books by J. Pedlosky |20j or R. K. Zeytounian [21]) there are few re- 
sults concerning the rigorous justification of the derivation that we presented in the previous 
subsection. To our knowledge, the first mathematical justification of Oberbeck-Boussinesq 
approximation in this context has been given only rather recently in the framework of the 
so-called variational weak solutions (see [H] for a complete presentation of such solutions for 
the full Navier-Stokes equations). The case of bounded domains with potential V = xs (or 
more generally, in W^''^{Q)) has been treated by E. Feireisl and A. Novotny in |12^ll3j. while 
the exterior domain case has been studied by E. Feireisl and M. Schonbek in [14] (still under 
the assumption V E W^'°°{Q,), thus ruling out the common but not so physical assumption 
that V = X3). For passing to the limit, all those works borrow some seminal ideas that have 
been introduced by P.-L. Lions in his book ^IQJ and B. Desjardins et al in flOJ in the related 
context of low Mach number limit for the isentropic Navier-Stokes equation^J- 

On the one hand, those results are very general for one may consider any finite energy 
data. On the other hand, the convergence results are not very accurate for they strongly rely 
on compactness methods : in particular convergence holds up to extraction only, and no rate 
may be given. 

1.3. Aim of the paper. Getting stronger results of convergence that is in particular con- 
vergence of the whole sequence with an explicit rate, is the main purpose of the present 
work. Considering general variational solutions is hopeless. We shall focus on strong solu- 
tions with the so-called critical regularity, a framework which is nowadays classical for the 
study of viscous compressible fluids (see e.g. [2j HI i5j). Of course, this will enforce us to 
restrict considerably the set of admissible data, but we will get much more accurate results 
of convergence. 

Working in a functional framework that has the same scaling invariance as (ll.2p . if any, 
is the basic idea. Here we see that (if T/^ = to simplify the presentation), the system is 
"almost" invariant for all ^ > by the rescaling 

a"(t, x) -^ a'{XH, Xx), u^t, x) -^ Xu'{XH, Ax), 9'{t, x) -^ X^e^{X\ Ax). 



For other recent results concerning the low Mach number asymptotics for the full Navier-Stokes equations, 
the reader may refer to [T] ITSl ITfl 118] . 
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If we believe in an energy type method then a good candidate for initial data is thus the 
homogeneous Sobolev space 

or rather the slightly smaller following homogeneous Besov space: 



which has nicer embedding properties {B2 i is embedded in bounded functions for instance) 
and better behaves with respect to maximal parabolic estimates. 

However, owing to the lower order pressure term, the above scaling invariance is not 
quite respected. Consequently, we have to work at a different level of regularity for the low 
frequencies of a^ and 9^, to compensate this scaling defect. All this is now well understood 
and already occurs in the isentropic case [3]. 

Finally, in the case k = that we shall also consider (and that cannot be studied in the 
framework of variational solutions), only the velocity is smoothed out during the evolution, 
and it is no longer possible to use a critical regularity framework: we will have to assume 
much more regularity. 

We end this introductory part with a short description of the rest of the paper. After an 
unavoidable introduction of some notations and functional spaces, the next section is devoted 
to the presentation of the main results of the paper. The analysis of the heat conducting case 
is carried out in Section [3] while k = is considered in Section HI Some technical estimates 
are postponed in the Appendix. 

2. Results 

Before presenting the main statements of the paper, we briefly introduce some notations 
and function spaces. We are given an homogeneous Littlewood-Paley decomposition {Aj)j^x 
that is a dyadic decomposition in the Fourier space for M. . One may for instance set Aj := 
(p{2~^ D) with (p{£,) := x(C/2) — x(C)i ^^^ X & non-increasing nonnegative smooth function 
supported in i?(0,4/3), and value 1 on i?(0, 3/4) (see |[2j. Chap. 2 for more details). 

We then define, for p G [1, +00] and s G M, the semi- norms 






j€Z 

In order to avoid complications due to polynomials, we adopt the following definition of 
homogeneous Besov spaces: 

S'l = |zG5'(M^) : ||z||r. <oo and lim SjZ = o} with Sj := xi^'^ D). 

To compensate the lack of strict scaling invariance of the system under consideration (as 
pointed out in the previous section), we also need to introduce the following hybrid Besov 
spaces with different regularity exponent in low and high frequencies: 

Definition 2.1. For s G M, p G [l,oo] and a > 0, we set 



•2 6s, ± 



Y^ 0.7s /„;„/„ -1 n.n\^l 






2J"(min(Q-\2J))^lAjz||LP 
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and define 

m't := (zeS' (M.^):\\z\\jis, ± < 00} and lim S'.z = oj- 

We shall mainly use the above definition with p = 2, in which case, the corresponding 
hybrid Besov space will be simply denoted by Ba , if the fact that p = 2 is clear from the 
context. 

We agree thatlj: 

(2.12) z^ := Y, ^J^ and z'' := ^ A^z. 

With this notation, we have 

\\Z\\ o^ii ^ h^ ds±1 -\- (X \\Z 73s 

" "-Dp, a " -"p,! -"p.l 

Therefore Bt, ^ is the bulk regularity of a function in Bp]oL while the behavior at infinity is 
given by the low frequency part which is in Bt^^ . Of course, changing the value of a does 
not affect the space, and the corresponding norms are equivalent. However a suitable choice 
of a will enable us to get uniform estimates with respect to e. 

As we shall work with time- dependent junctions with values in Besov spaces, we introduce 
the norms: 

""""-f^TC-B;,!) '^ IIII^'^*'")II-b;,iIIl'j(o,t) a^^*^ II^IIl?,(Bp;J') •" IIII^(*'")IIb;:±mIl'j(o,t)- 

As in many works using parabolic estimates in Besov spaces, it is somehow natural to take 
the time-Lebesgue norm before performing the summation for computing the Besov norm. 
This motivates us to introduce the following quantities: 

W'^WlUb;,) := J^2^1I4-^^IIl«.(L^) and ||n||^^(^.,±) := ^ 2^^(min(a-\ 2^))^^||AjM||i^(iP). 

The index T will be omitted if T = +00 and we shall denote by C{Bpi) (resp. C{Bp]a)) 
the subset of L°°{B^^) (resp. L°°{Bp]a)) constituted by continuous functions over R"*" with 
values in B^^ (resp. Bp\a)- 

Let us emphasize that, owing to Minkowski inequality, we have 

with equality if and only if g = 1. Similar properties hold for hybrid Besov spaces. 
Throughout, we shall denote 

(2.13) K := k/z^, a := X/v, Ji := ^/u with z^ := A + 2/i. 

One can state our first main result : the global existence of solutions corresponding to 
small (critical) data with estimates independent of e in the case k > 0. 

Theorem 2.1. Assume that the initial data (6q, Uq,9q) and that the potential term V^ satisfy, 
for a STnall enough constant i] depending only on k and fi: 

(2.14) 11^811 3,. + ||n§||. +\m--i^.<W, 

(2.15) i/5||vy^|| 3_ +\\V^\\ 3_ +\\dtV^\\ 3_ <7]u. 

L^m ) i-(B?; ) L^{Bj^ ) 



We omit the dependency with respect to the threshold a in the above notation because the value of a will 
be always clear from the context. 
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Let Oq := 6q + ^^(0). Then System ()1.2p with initial data (1 + eaQ,UQ, 1 + sOq) has a unique 
global solution (a^,u^,9^) (with a^ = ¥ + V^) which satisfies 

and, for a constant K depending only on k and /i, 

116^11 3_ +z^||6^|| 3 , + llu^ll 1 +i^|k^|| 5 +1161^11 _i , 

+z^lir|| 3, <K(\m\ 3 +\\u%\\ 1 +||6'g|| i,+\\dtV'\\ 3 ). 

Remark 2.1. Smoother data give rise to smoother solutions. For example if in addition to 
the above hypotheses, we have 

e\m\ 5 +u-^\\(eLu%)\\ 3 +e\\dtV''\\ s +e||Vy^|| s <r], 

then the above solution also satisfies 

eWb'W 5 +zy~^||(r,n^)|| 3 +eiy\\b''\\ s , + £i/||(u%r)|| r < Kt]. 

Next, combining this result with Strichartz estimates, we shah prove the following result 
of convergence to the Boussinesq system. 

Theorem 2.2. Consider a family of data (6Q,nQ,0Q, y^)e>o satisfying the conditions of The- 
orem \2.1\ with in addition 

Mo:=sup(||6§|| 3,_ + K|| 1 +110^11 1,, 

(2.16) ^>° '" , ^'^ 

+ z^2||Vy^|| 3_+||y^|| s_+\\dtV^\\ 3_)<r]u. 

L^{Bj. ) i°°(B?. ) LHB^: Y 

Let q^ := {9^ + W)l\f2 and Q" := {9" - b^)/V2. Assume that (Vuf^, Gg, V) converges (in the 
sense of distributions) to some triplet {vq,Qo,V) such that 

voeBl^, eoG5|i, VVGL^Bl^), dtVeL^Bl,). 

Then the following properties hold true : 

(1) System (jl.2p with initial data (1 + eOq^Uq, 1 + £9q) has a unique global solution with 
the properties described in Theorem \2.1\ 

_ . 1 .5 

(2) Boussinesq system ()1.7p admits a unique global solution {v,Q) in C{B2i) H L^{B2i) 
satisfying for some constant K = K{k, ft): 

\\{v,e)\\_ , +u\\{v,Q)\\ 5 <K(||(t;o,eo)||.i +\\dtV\\ 1 ). 

(3) The functions (f and Qu^ go to in the following meaning for all p € [2, oo] and 

s£[-l/2 + 4/p,3/p]: 






''^h'\\mB;-l-+) + ^^WQ^'WlHb;^,) < ^M'' Mo. 

(4) The couple {Vu^,Q^) tends to {v,Q) in the following meaning for all p € [2,oo] and 
s £ [-1/2 + 4/p, 3/p] with s > 1/2 : 

with dd' := e^ - e, dv' := Vu' - V, dV^ ■=V^ -V and C = C{fi, k, s,p). 
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Remark 2.2. If the data are smoother, e.g. as in Remark \2. 1\ then the results of convergence 
hold for stronger norms. For instance, it may be shown that {Qu^ ^q^") -^ in L^{Bp^ ^), 

4,1 4_3 i„l 4_3 

that Vu^ ^v in L^{B^^^ ') n L°°{B^^i '), and that Q^ ^ Q in L^{B^-^ ') n L°°{B^^-^ '), with 
the decay rate e^ p . 

Let us finally state our main global existence and convergence result for nonconducting 
fluids. 

Theorem 2.3. Assume that the initial data (ao,'Uo)'^o) ^'^^ ^^^ force term V^ verify that 
Cg:=||(aS,n§,7^§)|| 1 +(e^/)3||(a§,7^S)|| r + (ez.)2||nS|| 5 

^2,1 °2,1 ^2,1 

\\V^\\ 1 +(ezv)3||y^|| ^ +i/||y^|| 5 + u(eu)'^\\V'\\ 9 < w 

where the constant rj is sufficiently small and depends only on jl. 

Then System (jl.9p admits a unique global solution (a^,u^,7^^) which satisfies 

a^ G C{Bl^r\Bl^), u^ G C{B\^r\B\-^r\L^ (B\^r\B\^, n^ G c{Bl^r\Bl^)r\L^{Bh'^r\Bh'^) 

and, for some constant K depending only on jl, 

i°°(B2^i) L^iBl,) L^{Bl,) L^{Bl,) 

+^^||7^^|| 3 , +^^(e^^)2||7^^|| 7, +z^|k^|| 5 + z^(ez^)2||n^|| 9 < KQ. 

LHB7J+) LHB^J+) ii(B|,) i^Ba^i) 

Suppose in addition that 0g — )• Qq, that Vuq — )• vq and that V^ ^ V with 

(2.17) Iboll 1 +||Vy|| 3 +IIG0II 1 <??w. 

Then the corresponding limit Boussinesq system (11. lip admits a unique global solution {Q,v) 
in C(i32^]^) X {C{B2i) r\L^{B2i))- Furthermore we have 

(2.18) 11(6,^)11 1 +/.||^|| 5 < K||(eo,^o)|| . 1 . 

L°°{Bl^) L^Bl^) Bl^ 

In addition, if Cq is bounded by some constant Cq when e goes to then {Qu^ ,Vf) goes to 
zero with the following rates of convergence for all p ^ [2, 00) : 

(2.19) ||(Qn^7^^)|| 2p 21 <KCQe'^-'v 

(2.20) J^^II(Qw",^')llz2(iJ;_j < KCo{eu)l'' if s G [-1/2 + 4/p, 3/p]. 

Finally, if Qq and Vuq are independent of e then for all p and s as above (with in addition 
s > 1/2;, and T > 0, 

G"-e^O m CT{B;f) and Vu'-v^O in CT{B'j^J^+B'^f)n{Ll{B'p^^)+L^j.{B;^^)), 

--S 

and the rate of convergence is ep 

The above statements deserve some comments: 



The reader may refer to Inequalities (|4.116p and (|4.117p for the general case. 
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(1) In this paper, for simplicity, we focussed on the physical dimension 3. However similar 
statements may be established in any dimension d>2. 

(2) In the case of large data, we expect, as for the isentropic Navier-Stokes equations stud- 
ied in [6] , the lifespan of the solutions to (jl.2p to tend to that of the limit Oberbeck- 
Boussinesq equations. Global existence for the limit equations should entail global 
existence for (jl.2p with small e, if k > 0. This is of particular interest in dimen- 
sion two, as the limit equations are globally well-posed for any data with the above 
smoothness. We reserve this study to future works. 

(3) We also reserve the case of other boundary conditions, in particular the periodic ones, 
to future works. We want to point out that the global existence statements (that is 
Theorem 12.11 as well as the first part of Theorem I2.3P extend to that case. At the 
same time, no dispersive inequalities are available, hence the approach for proving 
convergence is expected to be completely different, provided based on the filtering 
method, as in the isentropic case [7]. 

We end this section by explaining the general strategy for the proof of convergence. The 
first step consists in proving uniform global a priori estimates. This in fact corresponds to 
the statement of Theorem 12.11 and to the first part of Theorem 12.31 We shall see that the 
proof reduces to the case e = 1 after suitable rescaling of the equations. Then, proving 
convergence requires two steps : first we establish that the oscillating part of the solution 
converges to (this relies on Strichartz estimates), and next establish strong convergence to 
Oberbeck-Boussinesq for the incompressible modes. Note that, owing to the fact that only 
small solutions are considered, we do not need to resort to bootstrap arguments. 

3. Global existence and convergence in the case k > 
Let us first notice that performing the change of unknowry: 

(3.21) (6, u, 9){t, x) := eib" ,u' ,e')ie'^i^t, evx) 
and the change of data 

(3.22) {hQ,UQ,eQ){x) := £{bl,ul,9l){evx) and V{t,x) := eV'ie^vt^evx) 

reduces the study to the case u = 1 and e = 1. Indeed it is obvious that {b^,u'^,6^) satisfies 
(jl.2p if and only if {b,u,9) satisfies the same system with e = 1, Lame coefficients (A,/i) := 
z^~^(A, fi) and heat conductivity k := i^'^k, provided the data have been changed according to 
(j3.22p . This change of variables has the desired effect on the norms that are used in Theorem 
12.11 For example, we have, up to a constant independent of e and z^, 

||6oll 3 = i^~"'"||fenll 3 ) ll'Woll 1 = ^^"'"ll'Wnll 1 ) l|6'oll 1 + = i^~"^ll^nll i +j 
II uii^^,,- II 011^^,^,- > II uii^^,^ II oll^^y II 011^-^.+ II Oll^-,.+ > 

\\VVit,-)Ls,. = e\\VV%e^ut,-)L^,. and \\dtVit,-)\\ _^ = s^dtV%s^,.t,-)L^^,., 

B{ Bei B{ B/v 

hence 



iVyll ,3_ =u-h\VV^\\ _3_ and \\dtV\\ ,3_ =u-^\\dtV^^ 



3 



■L2(Bf ) L^B^J ) L^B^' ) Li(Bl,' ) 

Consequently, in order to prove Theorem 12. 1^ it is suffices to consider the case z^ = 1 and 
e = 1. We shall resume to the original variables only at the end of this section, for getting 
the convergence results of Theorem [ 



Recall that u = X + 2fi 
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3.1. The linearized system. In the case e 
read 



1, the hnearized equations about (0,0,0) 



(3.23) 



r (9t6 + div'u = 0, 
dtu - flAu - (A + /i)Vdiv'u + \7{b + 9) = 0, 
dtO + divu- kA9 = 0. 



We aim at proving energy type estimates for (J},u,9). Roughly speaking, we shall exhibit a 
low frequency parabolic type smoothing for all the components of the solution whereas, in 
high frequency, only (n, 9) will experience a parabolic smoothing. As for 6, it will be damped 
with no gain of regularity whatsoever. Throughout our proof (which will require several 
steps) we shall also pinpoint where one has to work in different level of regularities to get the 
aforementioned features of the system. 

Let us first notice that the gradient terms in the velocity equation involve only the potential 
part of the velocity. More precisely, setting d := A~^divu (with A* := |-D|*) and w := Vu = 
u + V(— A~^)div n, we get 



(3.24) 



r dth + Ad = 0, 
dtd - Ad - A{b + 9) = 0, 
dt9 + Ad- RA9 = 0, 

^ dtw — jlAw = 0. 



As the last equation is the standard heat equation with constant diffusion, we focus on the 
proof of estimates for the first three equations. After localization by means of the homoge- 
neous Littlewood-Paley decomposition {Aj)j^i, the obtained system reads 



(3.25) 



{ dtbj + Adj = 0, 
dtdj - Adj - A{bj + 9j) = 0, 



dt9j + Adj - kA9j = 
with bj := Ajb, dj := Ajd and 9j := Aj9. 

Step 1: Basic Energy Estimate for (6, d, 9). Owing to the antisymmetric structure of the first 
order terms in (I3.25p . we readily get 



(3.26) 



\> 



2 

'illL2 



\d 



■JUL^ 



+ \\9,\\l,]+\\Adj\\l,+K\\A9,\\l,=0. 



Step 2: Improved Energy Estimate for (b, d, 9) . We want to track the decay properties of b. 
For that we notice that the auxiliary function Ab — d satisfies: 

dt[Abj - dj] + A{bj + 9j) = 0. 

Hence taking the L^ inner product with Abj — dj yields 

^l\\Ab, - d,\\l, + ||A6,-||i. + {A9j\Ab,)^, - {{b,+9,)\Ad,)^, = 0, 

from which we deduce that 
1 d 



(3.27) 



2dt 



Abj - d 



■J\\L2 



+ \\bj\ 



L2 



+ 



'illL2 



+ \\Abj\\l2 + {A9j\Abj)^2 + K\\A9j\\l2 = 0. 
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Putting (j3.26p and (|3.27p together, we thus get for any a > 0, 



(3.28) 



2di 



a\\d,\\l, + \\Ab, - d,]\\l, + (1 + a)||6,||i. + (1 + a)\\9,\\l, 

+ \\Abj\\l2 + {AOjlAbj)^, + k(1 + a)\\Aej\\l2 + a||A(ij|||2 = 0. 



|2 

Il2- 



Let us denote 

(3.29) ff := a\\dj\\l, + (1 + a)\\bj\\l, + ||A6,- - djWl^ + (1 + a)\\ej\\l„ 

(3.30) F| := ^\\Ab,\\l,+a\\Adj\\l, + (^k{l + a) - 0||Ae,| 
Then combining (j3.28p with the following Young inequality: 

5l|A(>,lli.+^ 

1 d 



(Ae,|A6,),^.|<^||A^,||i, + ^||A6,-||i„ 



implies that 

(3.31) 

Let us notice that 

Therefore, because 



2dt 



f^ + H^<0. 



ff = {a + l)\\{b,,dj,e,)\\l, + \\Ab,\\l,-2{Ab,\d,)L2. 



2\{Abj\dj)L2\<'^\\Abj\\l2+^\\dj\\l2, 



we have 

(3.32) (a-l)||d,||i, + i||A6,||i. </|-(a + l)||(6„0,)lli2<(a + ^)||d,||i. + ^||A6,-||i.. 

Let us first assume that k < 1. Then we take a = 2/k — 1 and (I3.32p thus implies that 



fl 



At the same time, we have 



k-mbj,d,,9,) 



--ii 



^J^"3)\\L2 



if k2^^ < 1, 
+ \\Abj\\l2 if k22j > L 



lL2 



2^m{bj,dj,ej)\\ 



L2 



if K22i < 1, 



k-^\\idj,9j)\\l2 + \\Abj\\l2 if k2^^>l. 
Therefore, one may easily conclude that for some (universal) constant c S (0, 1], 



bj , dj , Oj 



m\\ 



L2 



< 



522^ 



-ct\ 



\\{bj,dj,9,){0)\\L2 if 22^K<1, 
iA6j,(ij,0j)(O)||i2 if 2'^^k> 1. 



(3.33) 

(3.34) \\ikAb„d„e,)mL2<e 

Let us now assume that k > 1. Then we take a = 1 so that following the above computations 
after replacing everywhere k by 1, it is easy to conclude that 

c22Jt| 



(3.35) 
(3.36) 



\\{b„d,,ej){t)h2 
\\iAb,,d„e,){t)\\ 



< 



''\\{b„d,,ej)m\L2 if i<o. 



L2 



< 



-ct 



\\{Abj,dj, 



j,U,j,Uj 



L2 if J > 0. 



Therefore, denoting k = min(l,K) and putting together ()3.33p . ()3.34p . (j3.35p and (j3.36p . we 
end up with 

--2^^*||(6„d„0,)(O)||,.2 if 22iK<l, 



(3.37) 



\\{b„d,,e,)it)h2 



< 



IA6. 



'j[.T^)\\L2 



k-mdj,9j){t)\\L.<e-'^'{\\Abjm\L. + \\i^-Hdj,9jm\\L2) if 2^^k>l. 
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For reasons that will appear more clearly in the following steps, it is suitable to work with 
one less derivative in the high frequency regime. Now from the second inequality of (j3.37p 
and Bernstein inequality, we get for 2^-^k > 1, 

(3.38) ||K6,(t)||i. + ||A-i(d,-,0,)(i)llL2<e-^*(ll^&.-(O)||L2 + ||A-^d„0,-)(O)||LO- 

Step 3: Parabolic smoothing for 9. We here aim at tracking the high-frequency parabolic 
smoothing for 6. For that, we rewrite the last two equations of ()3.25p as follows 

J dtA-^dj - A{A-'^dj) - Oj = bj, 
[ dtA-^Oj - kA{A-^ej) + dj = 0. 

Then applying a direct energy method, we readily get 



1 d 
2di 



\\A-^d,\\l, + \\A-%\\1,^ + \\d,\\l,+k\\ej\\l, = ibj\A-~'d,). 

Therefore, performing a time integration yields 

\\A-\dj,ej){t)\\L2+ck [ \\A{dj,ej)\\L2dr<\\A-\d„9j){0)\\L2+ [ Wbjh^dr, 

Jo Jo 

and taking advantage of the second inequality of p.37p eventually leads to 

(3.39) \\A-Hdj,9,){t)h2+k [ \\A{d,,ej)\\L2dT<\\b,{0)h2+k-'\\A-\dj,ej){0)\\L2 

Jo 



in the high frequency regime, that is whenever 2^\/k > 1. 

Step 4^ Parabolic smoothing for d. Given that 

dtdj - Adj = A{bj + 0j), 
one may write that 

\\dj{t)\\L2+c2^^ I \\dj\\L2dT<\\d,{mL'2+ I \\A{b,,e,)\\L2dT. 

Jo Jo 

The previous steps ensure that, for 2^\/k > 1, 

' ||A6,-||^2dr < \\Abj{0)h2+k-^\\{dj,ejmh2, 
\\Aej\\L2dT < k-^b,{0)\\L2+k-^A~Hd„9,)m\L2. 



'0 

Therefore we have 



(3.40) 2^^ [ \\dj\\L2dT<\\Abj{0)\\L2+k-^\\bjm\L2 
Jo 

+ k-^\\A-Hd„9,)m\L2 + k-'\\{d„9,){0)U2. 
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Step 5: Final a priori estimate for {b,u,6). Putting together inequalities (|3.37p . (j3.39p and 
(j3.40p and using the standard properties of the heat equation (as regards w), we get if j < 0: 

(3.41) \\{b,,u„9,){t)h2+2^^ [ \\{bj,uj,ej)\\L2dT < C\\{bj,uj,9j)m\L2, 

Jo 

and, if j > 0: 

(3.42) \\{2^bj,Uj,2-^e,)it)\\L2+ [ \\i2^bj,2^^Uj,2^ej)\\L2dT<cm^bj,Uj,2-^e,mh2. 

Jo 

The above constant C depends only on fi and k. 

3.2. A priori estimates for the paralinearized system. As pointed out in the previous 
subsection (see in particular (|3.42p ). there is no gain of regularity for b throughout the 
evolution (only damping in fact). Therefore, the convection term v ■ V6 cannot just be 
considered as a source term, tractable by Duhamel formula, for the presence of V6 will 
induce a loss of one derivative in the estimates. 

At the same time, at the level of L^ estimates, this convection term is rather harmless 
provided divu is in L^(M^;L°°) (it is only a matter of integrating bv parts). The natural 
idea is thus to keep the convection terms in the linearized equationq^ and to resume to the 
method of the previous paragraph. As however the Littlewood-Paley localization operator Aj 
does not commute with the material derivative (9t + w • V), it is convenient to keep only the 
'bad' part of the convection term, that is the one which does induce a loss of one derivative. 
In order to better explain what we mean, we have to give a short presentation of Bony's 
decomposition (first introduced in [3j) and paraproduct calculus. The paraproduct is the 
bilinear operator defined on the set of couples of tempered distributions, by 

Tfg:=Y,Sj^if^j9 with Sj.i := xi2-^'''^D). 
j 

The (formal) Bony decomposition of the product fg reads 

fg = Tjg + T'gf. 

The basic idea is that the term Tfg is always defined but cannot be more regular than g, 
and that under suitable assumptions the other term T' f is more regular. If we look at 
the convection term, the 'bad' part that may cause a loss of one derivative and has to be 
included in the linear analysis is thus (with the summation convention over repeated indices) 
T^kdkb. This motivates us to extend the analysis of the previous subsection to the following 
'paralinearized' system: 



(3.43) 



dtb + M + T^kdkb = B, 
did + T^kdkd -Ad- A(6 + e) = D, 
dtO + M + T^kdkO - RAO = G, 
^ dfW + T^kdkW — flAw = W, 

where the source terms B, D, G, W and the vector field v are given. 



We keep all the terms just for questions of symmetry, but only v ■ V6 may cause a loss of derivative. 
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Proposition 3.1. Let V(t) := /q ||Vf ||loo dr. For all s G M, there exists a constant K 
depending only on fl, k, and a universal constant C such that the following inequality holds 
true: 

t V 1 / t \ "^il' t V 1 / / 1". 1 z.l 1 






Proof. Compared to the study of the previous subsection, the main additional difficuhy lies 
in the paraconvection terms. Indeed, the source terms may be easily dealt with by means of 
the Duhamel formula. 

The paraconvection terms may be handled thanks to the following inequality: 

(3.44) \{ct>{2-W){T,udkz)\(t){2-W)z)^,\<C\\Vv\\L^\\(t>{2-W)z\\L2 ^ U{2-^' D)z\\l2 

\r-j\<N 

which holds true for any smooth function (f) with compact support away from the origin and 
large enough integer A^ depending only on Supp cj) and Supp ip. 

Let us justify ()3.44p . We fix some integer A^ so that 

Supp0(2~^-) n Supp (x(2"^"-) * V'(2"^"-)) = whenever \j - j'\ > N. 
Then we use the following algebraic identity: 
{<l){2'iD){T,,dkz)\<l){2'W)z)^, = Yl {(l>{'^''D){Sy-iv''dkAfzM2-W)z)^, 

\J'-J\<N 

= Yl {H'^''miSj'-i-Sj-iydkA,,zM2-W)z)^, 

\J'-J\<N 

+ Yl {m2~W),Sj-,v'']dkAj,z\(l)i2-W)z)^, 

\i'~j\<N 

+{Sj.iv^dkcl){2-W)z\<P{2-W)z)L2. 

The first term may be bounded thanks to spectral localization and Bernstein inequality, and 
the second, to a standard commutator estimate (see e.g. [2], Lemma 2.97). The last term 
may be dealt with according to the following integration by parts: 

I Sj-iv^dk<t){2-^ D)z <t){2-^ D)z dx = -]- f divSj-iv {<i){2~W)zf dx. 

Let us now resume to the proof of Proposition 13.11 As an example, we show how the first 
two steps of the previous subsection have to be adapted for ()3.43p . So we apply Aj to the 
first three equations and get: 

dthj + Mj + /^j{T.^kdkh) = Bj, 
dtdj + Aj{T^kdkd) - Adj - A{bj + Oj) = Dj, 
^ dtOj + Adj + AjiT^kdkO) - kA9j = Gj. 
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Taking the L^-inner product of the first, second and third equations with bj, dj and 9j, 
respectively, we find that 

Ijt^Mh + \\dA\h + ll^.lli^] + IIAdilli^ +^||A^,-||i. + {AjiT^^d,b)\Ajb)^, 
Therefore using Inequahty ()3.44p we readily get 

~m,d„e,)\\l + \\Ad,\\l, + k\\A6,\\l, < \\{b„dj,ej)h2 

x(\\{B,,Dj,Gj)\\l2+C\\Vv\\l^ Y, \\ibj',dj.,ej.)\\L2). 

\f-J\<N 

Next, we use the fact that Abj — dj satisfies 

dt{Abj -dj) + A{bj +dj) + AAj{T^kdkb) - Aj{T^kdkd) = ABj - Dj. 
Therefore arguing as in the second step of the previous section, we get 

—ff + H] + ((AA,(T,.9fc6) - A,{T,.dkdmAb, - d,))^, 

Hl + a){^j{K,dkb)\Ajb)^,+a{Aj{T^kdkd)\A,d)^, + {l + a){A,{T^kdke)\A,9)^, 

= (1 + a){Bj\bj)L2 + a{Dj\dj)L2 + (1 + a)(G,|0,)i2 + {{AB^ - Dj)\{Ab, - dj)) ^, 

where f^ and H^ have been defined in (j3.29p and (j3.30p . and a = 2/k — 1. 

Note that all the paraconvection terms except the first one may be directly dealt with 
according to (j3.44p . As for the first one, we may use the decomposition: 

AAj{T,udkb) - Aj{T,udkd) = A,T,.9fc(A6 - d) + 2^' [0(2-^1)), T,u]dkb 

with (j){S,) := |^|v(0- Therefore, applying again (j3.44p and Lemma 2.97 in ^, we end up with 

\{{AA,{T,kdkb) - Aj{T,,dkd))\{Ab, - dj))^,\ 

< \\Vv\\L^\\Abj -dj\\L2 ^ {\\Abj' -dj>\\L2 + ||A6jv||^2). 

The following steps may be done similarly, once noticed that operators such as A Aj 
may be written 2^^ (j)(2~^ D) for some suitable function <p with the same support as y?. The 
final inequality may be obtained after multiplying by 2^^, performing a summation over j 
and applying Gronwall's lemma. The details are left to the reader. D 

3.3. The proof of global existence. This paragraph is devoted to proving Theorem 12.11 

in the case e = z/ = 1. As explained at the incipit of this section, this will imply the global 

existence for general positive e and z/. The proof of existence and uniqueness is similar to that 

for the full Navier-Stokes system in [5j. The only difference here is that the source term W^ 
_ 1 _ 1 

is not in L^(W^; B2i)- However it still belongs to Lj^^{W^ ; B2 i) which suffices to establish 
local- in-time results, global results being a consequence of the following a priori estimates. 
Note that a direct proof based on Priedrichs spectral truncation method may also be easily 
implemented as we are interested in L^ type estimates. 
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So let us now derive global a priori estimates under the smallness assumptions (j2.14p and 
(j2.15p . Such estimates rely on Proposition 13.11 with s = 1/2, once noticed that 

u = T>u + (Id -V)u = w- VA^^d, 



that (6, d,9,w) satisfies ()3.43p with v = u and, using the summation convention over repeated 
indices, 



B ■= T^kdkb-u-Vb-bdivu-dtV-div{uV), 
D := T^kdkd-A'^diY{u-Vu)-A-^dw 



a /- » /T ~N^,. N (0 — a)Va 
{fiAu + {X + fi)Vdivu)+ ^ ' 



1 + a 



(1 + a) 



G ■= TkOkO-u-ve-edwu 



1 



-RAO + ^^[2/2|Z?mP + Afdivu)^], 
1 + a 1 + a 

« /- * /T ~^^ ,. N (^ ~ a)Va 
3^(M- + (A + /.)Vd.vn) + i^^ 



W := T^kdkW-V{u-Vu)-V 
Setting {/(t) := /^ ||Vn||Loo dr and 

X{t):=\\b\\_ 3_ +||mL 1 +11^11- - 1+ + / (ll^ll 3 + + ||7/|| 5 +11^11 3^)dT, 



we may write 



(3.45) X{t) < Ke^^W('x(0) + / e-^^M(||S||.3,_ 



{D^W)\\^.^^ + \\G\\^_.^,)dry 



Throughout we suppose that 1 + a is bounded and bounded away from 0, an assumption that 
is satisfied provided 11 a 11 3 is small enough. 

L^iBl,) 
Bounding \\B\\ _i,3 . According to Bony's decomposition, we have 

u-Vb-T^,dkb = T'Q^,u\ 

Hence standard results for the paraproduct imply (just decompose b into low and high fre- 
quencies) : 



(3.46) 






BT 



Likewise, according to Lemma |5. 11 we have 



(3.47) 



|6divM|| 3 < ||6|| 3 lldiviil 
I 11^^^, _ _ II ii^^,,_ii I 



3 . 



Finally, because div (yu) = T^diva + u ■ W, we have 



(3.48) 



div (uy) 3 < U^ 3 divu 3 + Vy 3 \\u\\ 3 . 

Bi i'2,1 ^1 B21 



■B?- 
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Bounding \\(D,W)\\ . i . We concentrate on D, proving estimates for W being similar. We 

^2,1 

have 

T^kdkd - A-Miv {u ■ Vu) = [T^k,A-^di]dku' - A-^diTi^^^.uK 

Therefore, resorting to standard commutator estimates and continuity results for the para- 
product (see e.g. [2]), we get 



(3.49) 



\T,,kdi.d — A div (n • Vn) i 5, Vn roo n i . 

\ U K \ n\ r^ \\ \\U II II 

^2,1 ^2,1 



Next, combining composition and product estimates yields 



(3.50) 
and also 



l + a 



-V^u\ 



< 



^2,1 



\a\\ . 3 \\u\\ 



R 2 R2 

■^^2,1 ^2,1 



a)Va 



l + a 



1 ^ (1 + ll«ll 3 )l|Va|| 1 I 

^2,1 -°2,1 ^2.1 



I .3 + 
-^2.1 



3 + \\a\\ 3 ). 

■t*2,l ^2,1 



.3 .3 

Note that we expect 9^ and 6^ to belong to L^ (M+ ; i?! -^ ) and L-'^(]R+; i?2^-^), respectively, and 



that, applying Holder inequality yields 
9-a)Va 



l + a 



< 1 + 1 



~ 3 _ .3 

So finally, because B^ ' "^ B2 i , 
a)Va 



a 3 



{a,e^)\\ 3 +||o|| 3 ii^^ii 3 ). 
'"l^{bi,) "L^(Biy L^{Bi,y 



(3.51) 



l + a 



Li{B|J 



< 1 + 



L'=°{B^' Y 



1 + + lla| 



Bounding \\G\ 

Therefore 

(3.52) 



We first use the fact that 



\T„kdkO-u-Ve\\ 1 , < ||V6i|| 3,||u|| 5. 

^1 ^1 ^2,1 



L°°{B^' ) Li{Bi^' ) 



0- 



Next, Lemma |5 . II implies that 
(3.53) 

(3.54) 



l^divul 



< 



Bi 



1 4_ dlVM 3 , 



l + a B^ 



< 



B2 R^ 



Finally, since B2 I 
(3.55) 



M+ 



-B^ ^' , standard product laws enable us to write that 



1 



l + a 



Vu0Vu 



1 + ^ (1 + l|a|| 3 )\\Vuf.u2- 
^-,.+ ~ V II 11^3^711 ii^va 



Plugging inequalities ()3.46p to p.55p in (I3.45P and making the assumption that 



(3.56) 



iVnl 



Li(L°°) 



< 1 and ||y|| 



L°°{Bf' ) 



+ livy| 



LHB^ ) 



<1, 
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we thus get 

Xit)<C(XiO) + \\dtV\\ 3_ +X\t) + X\t)). 

It is now clear that the solution may be bounded for all time if ^(0) and V are small enough: 
we get for some constant K depending only on k, fl and A, 

(3.57) X{t) < KCo 

with 

Co := ||6o|| 3.- + ll^oll . 1 + ll^oll --1,+ + ll^iV^II ^ 3^_ . 

B^' Bl^ B^ ^ L^{B^ ) 

3.4. Convergence to the viscous and diffusive Boussinesq system. The key observa- 
tion is that in the asymptotics e going to 0, the leading order part of the system for (q^ , Qu^) is 
the acoustic wave equation, which has dispersive properties. This will enable us to show (first 

step) that {q^ , Qu^) tends strongly to in some negative Besov space. Next, we shall check 

. 1 . _i 

that the limit Boussinesq system (ll.7p supplemented with small data vq £ B2 ^ , ©o S i?2 1 

.1 . j_ 

and potential V with dtV G L^{B2i) and W G Lp'{B2i) has a unique global solution. Fi- 
nally, resorting to maximal regularity estimates for the heat equation, we will conclude that 
(PM^G^)^ {v,Q). 

3.4.1. Convergence to zero for the oscillating modes {q^, Qu^). In order to exhibit the decay 
properties of ((7^, Qu^), we only have to consider the case £ = 1 and ly = 1 thanks to the 
rescaling (j3.2ip . which implies in particular that 

{q, Qu){t, x) = e(g^, Qu^){e'^ut, evx). 

Then using Strichartz estimates for the acoustic wave equation (see Proposition 15.11 in the 
appendix) will enable us to bound some suitable norm of [q, Qu). Resuming to the original 
variables, we then get for free the convergence to for (q^ , Qu'^), with an explicit rate. 

Let us give more details : (g, Qu) satisfies 

dtq + V2div Qu = -div (qu) -—(dtV + div (Vu) + k ) 

2 V 1 + ay 



(3.58) 



+ ^j^[2f,\Duf + Xid\vuf] 



dtQu + V2Vq= qI 1'^ )Va--^-n-VnV 



Therefore Strichartz estimates (first inequality of Proposition 15.11 with s = 1/2) enable us to 

, 2p .2-1 

bound the norm of {q, Qu) in Lp-^ i^pi ^ ) for all p G [2, 00) in terms of the norm of the data 

.1 . j_ 

in B21 and of the right-hand side in L^{B2i)- Under our present assumptions however, the 



1 



last term in the r.h.s. of the first equation belongs only to the larger space L^{B-^ ^' ). So 
one has to use the second inequality of Proposition 15.11 and just get estimates in the wider 

space L^(S|;^^'^). 

__i _|_ 
Let us bound the r.h.s. of (j3.58p in L^{B^ ^' ). All the terms may be dealt with by taking 

advantage of standard product laws and Lemma 15. li More precisely we have, keeping in 
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mind the smallness of a in L°°{Bl' ) (and thus also in L°°{B^^) and L°°(M+ x M^)): 






< 



LHBJ) 
L 



{l + a)-^Au\\ 
u ■ Vu 



LHBl,) 

< 



V\\ 3 ||u|| 3 , 

L2(B|,) L^iBl,) 

^f 3 , 

L^Bl,) 



+11^ 



.^A.+J + ll«l 



^11 3 • 

LHBl,) 



L^Bl,) L^Bl,) L\Bl^) L--{Bl ) L^B^^) 

L'{BlS 

1 < Hull 

L^'.Bl^) 

Given that Q is a 0-th order multipher (hence maps all Besov spaces involved here into 
themselves) , that B2 l and i?! 1 ^'^s continuously embedded in i3^ ^ ' , and that B^ ' is 

continuously embedded in i3^ ^' , we eventually conclude that (with the notation of (|3.45p 
and (13371) ): 



\\{q,Qu)\\_ 2p 2_3 



iP^(BP^2 ■) 



, < lite. C..o)ll,- . . + X + A- + ll«.''ll,,,|.-, + ll^-ll,,,,!/- 



Therefore, given that X{t) < KCq and that Co is small, 



(3.59) 



\\iq,Qu)\\ 2p 2_3+ < KCo for aU p G [2, 00) 



with ET depending only on p, k and jl. 

On the other hand, Inequality (I3.57P implies that 



{q,Qu)\\ . <KCo. 
Li (B/ ) 



Therefore, using the fact that 



we get also 



[L\Br^),L—-{Bl^'' )Up+2)C L\Bl,') with g = (p + 2)/2, 



q,Qu)\\^ 4_3, <KCo for all gG[2,oo). 
L2(B« ^' ) 



Given that (g, Qu) is in L [B^' ) hence in L {B ^ "^ ), an ultimate interpolation ensures 



that 



(3.60) 



{qMu)\\i^^^^s-i.+^<KCQ forah s E [-1/2 + 4/p, 3/p] and p G [2,oo). 



Of course, we also have Qu in L'^{B2i) whence in L'^{B^^) for all p > 2. Therefore, interpo- 
lating with ()3.60p . we deduce that 



(3.61) 



\Qu\\i2(B^ ^ < KCo for ah s G [-1/2 + 4/p,3/p]. 
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Now coming back to the initial variables, (j3.59p . (j3.60p and p.6ip translate into 

(3.62) iy^'h\{Q',Qu')\\ 2P 2_3^ <K{eu)^-pC^ for all p G [2,oo), 

(3.63) u^\\q''\\^2^^s-i.+^<K{eu)p''C^ for all s £ [-1/2 + A/p,3/p], 

(3.64) u^\\Qu''\\i2(^B.^)<K{eij)l''C^ for all s G [-1/2 + 4/p, 3/p]. 

3.4.2. Global existence for the Boussinesq system (jl.7p . Let us first briefly justify that, under 
our assumptions, the limit data (Go, vq, V) give rise to a global solution to ()1.7p . Establishing 
this is an obvious modification of the proof for the standard incompressible Navier-Stokes 
equation. It is only a matter of rewriting the system as 

e(i) = e*t^eo+ f e^^-^'>i^(^{dtV + v-VV)-v-Ve]dT, 

v{t) = e'^^vo - I e^(*-")^p(7; • Vv + ^©Vf) dr, 

and the global-in-time solvability for small data may be achieved as a consequence of the 
Banach fixed point theorem. Let us just check that global a priori estimates are available in 
the case of small data. Applying Proposition 15.21 and using that the product is continuous 

.1 .3 .1 

from B21 X B2^ to B21 implies that 

lieii 1 +K||e|| 5 < iieoll i +\\dtV\\ i 

+ \\v\\ 3 (IIVGII 1 +||Vl^|| 1 ) 



and that 



I 1 +/i||t'|| 5 5: llfoll 1 + lli'll 3 ||Vf|| 1 + ||Vy|| 1 ||0|| 3 . 

L^iBl,) LHBi,) B|, L2{Bl,) L^Bi,) L^Bi,) L^Bi,) 



Hence, setting 



Y:=\\{e,v)\\ ^1+^11(6,^)11 5 



we get for some constant K = K(fi,K) 



Y <K(Yo + \\dtV\\ 1 +Y(Y + u-2\\VV\\ 1 )), 



1 



and it thus easy to close the estimates globally if In, 119/^11 1 and u^ IIVGII 1 are 
small compared to i'. 

3.4.3. Convergence for the "incompressible" modes {Q^,Vu^). In this paragraph, we prove 
the convergence of (0^,P«^) to the solution (0,w) to the Boussinesq equation (jl.7p . We 
claim that for any p G [2, 00) and s G [—1/2 + 4/p, 3/p] with s > 1/2 : 

• gQs ._Qs _Q tgn^ig tQ Q ij^ L'^(Bp~y) n L°^{Bp~y), 

• dv^ ■= Pu^ - V tends to in L^Bpltu) n L'^iBp'ey). 
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For proving that, we shall use the parabolic estimates of Proposition 15.21 for the system 
satisfied by {dQ^,dv'^). Let us first focus on (59^. By performing the difference between (jl.5p 
and (jl.7p . we see that 

dtdd' - -Add^ = -Vu^ ■ V(5B" - *^ • Ve + —{dtS^^ + Vu^-VS^^ + dv^ ■ VV) 

+div {{V^ - e^)Qun + f Ag^ - ^K^-i^Ar + ^-^[2^|Z?^^|2 + A(divu^)2]. 

Hence, according to Proposition [521 it suffices to get suitable estimates for the right-hand side 
in L^{Bp^eu ) + Lp'{Bp^eu )■ From product estimates (see Lemma [5?T]l we easily get under the 
assumption that s > 1/2 (in fact here we just need s > —1/2 owing to div dv"^ = divVu^ = 0): 

(3.66) 11*^- ve|Li,Bs-2.+, < iiveii 1 11*^11,2, B.-l.+^, 

(3.67) ll^^i'-V<5V^"||M.fi.-2,+. < llPn^ll 3 \\VSi^'\\j^,f.s-2,+., 

^ ' II ll-t^ i-Dp.ei^ ) ~ II L2(^2 ■) " "-t^ (.-Dp.ei/ }' 

(3.68) 11*^- Vy||M,As-2.+ , < ||Vy|| 1 ||&^|L2,5.-1,+ ,. 

We split the next term into (referring to the notation introduced in (j2.12p with a = ev^ 

div((y" - e")Qn") = div((y^ - G^'^)Q'u^) - div(e"'''Qu^). 
First we have 

||div((y^ - e^'^)Qn^)||,i(^.-..) < W^ - e^'')Q-^ll^.(B--) 

(3.69) <{W% 3 +||e^'^|| 3 )||Qn^|U,^.-i,+,, 

and, second 

■(e^''^Q7x^)iL,,B.-2,+. <— lle^•'^Qn^l^,.o»-l^<llQu"l^2.6. .\w^\ _^ 



(3.70) l|div(e-'^Q.^)||,,(^.-2,.) < ^lie-'^Qni,,^^.-.) < II2-1l2(b;,,)II0^''II,,^^^,.^ 

Next, we see that, for all a G [0, 1), 



I Aril 1 „, <|ka^|| 3 JIAril 1, . 



Now, by interpolation 



II --J 

-a rs-/ 



|a II . 3 - iiu, II 3 



< ll/i^l|l~°lln<^ll" 

ri; " 3 " . 1 



■ -^ — OL rs_/ o . i 

^2^1 ^2^1 ^2^1 



. 3 



and the definition of the norm in i?|J implies that 



la'^ll 1 + ez^lla^ll 3 < ||a'^|| 3 

I llo^ II 11^5 ~ II llo2' 

^2,1 ■°2.1 ^2,e 



Therefore 



(3.71) We^a^W 3 , < (ei^)"||a1| 3 . 

\ ) II II pj-a ~ V ; II II -5,- 

■°2,1 ■°2,£i/ 

~_i_Q,^4. ~-^-2-a,+ 

We also notice that B2 ly ' M- Bp^^u for p > 2. Therefore if we take 

a := 3/p — s, 
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then we get, keeping in mind that ||a^|| 3 _ is small, 

(3.72) | L ^°' M'\\j,,f,.-2.+,<v-\euy\\a'\\ 3 ||6'"|| 3,. 

Finally, 



e 



■[2/x|Z)n^|2 + A(divn^)2]|| ^^_, <e(l + ||ea^|| , a JUVtx^n^ 



1 + ^a^'" ' "■ ' '"l^bJ)^"'' ' "" "L--(Bl,y"' "l^bIS 



1||„£|| . Ml„,e||2 



<e{l + u-'\\a'\\^ ,^,.)\\u' 



3 



i°°(S2%) L^B^,) 



At this point, let us notice that for all z £ B21 and a G [0, 1], 

K ■ -i-a,+ ^ ■ i-a I l.t'^J K ■ - 

■"2,£i. ^2.1 ^2 



< ie,.r-'\ 



z\ 



' - ' ~,s-2, 



Since i?2 ^i, ' ^^ Bp^ev (with a = 3/p — s), we thus end up with 

(3.73) ||^^[2;.|Z)n^|2 + A(divu^)2]||^^^._2,+ <^-i(e^)"||n^f 3 • 

So putting ()3.65p to (j3.73p together and using (j2.16p . we conclude that 

^^ll<»'llL2(Bri'+) + \\^"h^iB;-i'+) ^ II«^oIIb-j.+ +Moll(*^«^)llli2(B;-i.+) 

(3.74) +Mo{eur{Mo + 1) + Mo||V5V-^ 11^,(^.-2.-.) + \\d,dV^\\^,^^s..,^^^^,^^.-s,+y 

Let us now concentrate on the proof of estimates for 6v^. We have, subtracting (|1.7p from 
(|1.5p and using (ll.6p . 



/2 
dtdu' - nAdv' +V{Vu^-Vdv' + &"-Vv) = -^P(e^V(5K^+(B"Vy+ q'VV - 2q'Vb') 

\ 1 + ea"^ 1 + ea^ 

Therefore, according to Proposition 15.21 and to the fact that P is a self-map on any homoge- 
neous Besov space, we have 

\\^mfoo(fis-2.+-. + i^||(5i;^|Li.fjs,+ ^1 < ||&;oIIr'>-2.+ + WPu"^ ■ V&;^Ll.R^-2,+^ 

+ 11*' • Vt.|||^,(^.-2,+ ) + ||e^V5K^||^i(^.-2,+ ) + ||c©^Vy||^i(^.-2,+ ) 

+ ||Qn^ • VPu^||^w^.-2,+^ + II-— -Au''\\^-,,^s~2,+. + II — — — Va^||^i.^s-2,+-^. 

The following inequalities stem from product laws (see Lemma l5.ip . under the assumption 
that s > -1/2: 

(3.75) ll^w'-V*^Li.R.-2,+. < llPn^ll 3 ||V&^|L2fBs-2,+., 

(3.76) ll'5y'^ • Vf III rl/fis-2,+ ^ < ||Vf II 1 ll&I^ILo/Bs-l.+ x- 
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Next we have, if s > 1/2, 



(3.77) 
(3.78) 
(3.79) 
(3.80) 
(3.81) 
(3.82) 



k'-VQu^lLi 



^ {^p,£iy 



s-2,+ 
V 



ll«)^vy|Lwp.-2,+ 

1 1 ^£ V7 ?,£^ 1 1 



< 



1 II Qn^" 



©11 1 + 
WW 1 ||(5B" 









VF^II 1 I 
V6^|| 1 ||( 

LHBh) 



li2(B^-i'+)' 



li2(B;-i'+)- 



So arguing as in the proof of (|3.72p . we get 



(3.83) 



Finally, 



ea 



l + ea= 



-Au^ 



< ;y-Vp;y^"lla^l 



^l,RS-2,+^ ^ 1/ [eh') \\a II _3 _ ii'u- II 5 



^ ^Va^ 1 „, < Va^ 1 



Hence using again that i?2 Jj ^-^ -Bp|^,y' and (|3.7ip . we conclude that 
(3.84) i.-^C'-^)' 



< z.-H6z^)"lla^l 



'l-(B,1' ) L^Bl,) L^Bi;l) 



,1 i + a 



LHBl,) 



So putting together inequalities p.75p to (j3.84p . we end up with 

^II*'IIli(b;;+ ) + \\^'h^{B;-l-+) ^ II*oIIb-2.+ 

+Mo(||(&;^«^)||^,(^.-l,+) + ||V(5V^^||^,(^.-i)) + (eurMlil + ^.^IMq). 

Bearing in mind p.65p . we thus see that if Mq is small enough with respect to v, 



z/2||(5B" 
(3.85) 



^l^(b;-1'+) ^ W'^' \\L^iB;zl'+) + '^W'^' \\l^(b;:+,) + II'»^'IIl^(b;-2.+) ^ ll('»o>o)b. 



-l,+ ^ + ||(©'^||?™,5s-2.+ N + 2^11*^11 r 1/5=,+ N + 11*^ 
t2 



2,+ 

'p,£U 



+M^{eiy)'' + Mo\\VSV' 



lHb;;,''+) + ll^*^^llLi{Br.''+)+i'(^ri'+) 



whenever s > 1/2, 4/p — 1/2 < s < 3/p and 2 < p < 00. This completes the proof of the 
theorem. 



3.5. The case of smoother data. In order to improve the results of convergence (see 
Remark 12. 2p , we need to have higher order a priori estimates for the linear system ()3.43p . 

, 2p .1-1 

In effect, if we want to have convergence in (j3.62p for the norm Lp-2^^p ^-j j-^^j^gj. ^j^an 
Lp-2 (Bp^su' ) then we need 6 to have the same regularity as b, namely B21. So we need in 



addition that ^o S B^' and, owing to linear coupling, this will enforce us to take uq G Bl 
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Here we just point out what has to be modified to our previous arguments so as to handle 
such data. Let us start with (j3.25p . We concentrate on the high frequency regime. First we 
notice that 

dtd - kM = -Ad. 

Hence standard energy estimates ensure that 

\\Ae,{t)h2 + K22i||A0,||^i(^.) < \\A9jiO)h2 + ||A2d,-|Li(^2). 
Taking advantage of (I3.42p . we thus get 

(3.86) 2^'||^.-(*)IIl2 +2=^i%||ii(i2) < C\\{2^bj,dj,2^ej)mL^- 

We also need more regularity for (6, d). This is given by (I3.42p after multiplying by 2^: 

(3.87) \\{2^^b„2^d„0,)mL^+ I \\{2^'b„2'^dj,2^^e,)\\L2dT<C\\i2^^bj,2^dj,0,)m\L^- 

Jo 

Arguing as in the proof of Proposition 13. H we thus deduce that 

||6||?oc.Bf+i--nsf+2'-) + \\id,w,0)\\i^,^s+i,-. + / (||6||^.+i,+^^s+2,+ + \\{d,w,e)\\^s+3,-) dr 

tvl 1 ' t \ i_ ' j r\ 11 1 



< 

ft 



I ii "ii^i nsf 



+ ^ e-^^W (||i?||^.+i,-n^.+2,- + II {D, W, G) 11^.+!,- ) dr) • 



Starting from this inequality and following the computations of Subsection 13.31 it is easy to 
get the result of Remark 12.11 Next, resorting to the first inequality of Proposition 15.11 with 
s = 1/2 and to nonlinear estimates, we get Remark | 



4. The nonconducting case 

As pointed out in the introduction, in the case k = 0, it is easier to work with TZ^ . 
The reason why is that the linearized equations for [u'^ ,1Z'^) are the same as those of the 
classical barotropic Navier-Stokes equations (see next paragraph). Apart from this purely 
technical point and the fact that one has to work with smoother data, the overall approach 
for investigating the global existence and low Mach number issues is the same : first we 
perform the change of variables 

(4.88) {a,u,n){t,x)=£{a\u'',n''){e^vt,evx) and V{t,x) = eV'ie^vt^evx), 

so as to reduce the proof of existence to the case e = v = 1, and next we take advantage of 
dispersive properties of the acoustic wave equation, and of parabolic estimates to establish 
the convergence to some suitable solution of the Boussinesq system with no heat conduction 
(namely pTTT]) ). 

4.1. Linear and paralinear estimates. If we decompose, as in the heat-conducting case, 
the velocity field u into its (reduced) potential part d, and its divergence- free part w, then 
the linearized system about reads 

dta + Ad = 0, 

. ajd - Ad - A7^ = 0, 

dtw — flAw = 0. 
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As in the heat-conducting case, w just fulfihs the heat equation. Next, we notice that (7^, d) 
satisfies the hnearized equation for the compressible modes of the barotropic Navier-Stokes 
equations. Hence, fohowing the method of [1], we gather that for some universal constant C, 

||(7^„d,)(t)||L2 + 22^• / \\{n„d,)\\L2dT<c\\{n„d,mh2 if i<o, 

\\{2m„d,)it)h2+ f \\i2^n„2''^d,)h2dT<c\\i2^n„d,m\\L^ if i>o. 

Jo 
Now, from the first and last equations of (|4.89p . we see that 

aj{t) - Ujit) = aj{0) - TZjiO) for all t G M+ . 

Hence, taking advantage of the above estimate for TZj, we get 

max(l,2^')||a,(t)||i2 < C(max(l, 2^■)||(a,•(0),7^,(0))||i2 + ||d,- (0)11^2). 

From those inequalities, arguing as in the case k > 0, one may deduce a priori estimates for 
the following paralinearized equations: 

dta + Ad + T^kdka = A, 
dtd + T^kdkd -Ad-An = D, 
^^■^°^ ^ dtU + Ad + T^kdkTZ = R, 

dtw + T^kdkW — jjAw = W, 

where the source terms A, D, R, W and the vector field v are given. 

More precisely, we have 

Proposition 4.1. Let V{t) := L \\Vv\\loo dr. There exists a constant K depending only on 
jl and a universal constant C such that for all s ^"K, the following inequality holds true: 

\\{a,n)\\l^^^^s+l.-■^ + \\{d,w)\\ioo^^Bi^-^ + / (||((i,u;)||^s+2 + \\n\\^s+i.+ )dT 



+ 1^ e-^'^^^WHAR^^s^.,- + \\{D,W)\\^JdTy 



4.2. The proof of global existence. Here, in the case e = i/ = 1, we want to prove the 
existence of a global solution (a, u, TZ) to (jl.Op with 



a G C{Bi,nBi,), u G C{Bi^r\Bi^)r\L'{Bi^r\Bi^), n G C{Bi,r\Bi^)r\L\Bi''r\Bi'' ). 

For that, this is mainly a matter of proving a priori estimates in this space, taking for granted 
the existence of a solution. Indeed, the a priori estimates that we are going to prove below 
would be the same for the system truncated by means of the Friedrichs method (see e.g. [2], 
Chap. 10 for the related case of the barotropic Navier-Stokes equation). 
More precisely, we have to bound: 

(4.91) X := ||(a,7^)|| i 7 +11^11 i 5 +11^11 s 9 + 1|7^|| 3, 7,. 

L°^(B^^nB^^) L^{B^^nB^^) L^iBl^nBl^) LHB^'^nB^' ) 
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As we have in mind to apply Proposition 14.11 (twice: once with s = 3/2 and once with 
s = 7/2), we rewrite (jl.Op as fohows: 

( dta + T^kdka + Ad = A, 

dtd + T^kdkd -Ad-An = D, 

dtU + Ad + T^kdkU = R, 
^ dtw + T^kdkW - iJ.Aw = W, 



(4.92) 



where 
A := T^kdka — u -Va — adivu, 

D := r„^afc(i- A"Miv(n- Vn) - A^Miv 



■(/iAu + (A + /i)VdivM) 



aV{n + V) 



1 + a ^'' "" ' " " ■ '' ' '' '' (1 + a) 

R ■= T^kdkn-u-Vn-TZdi-vu-dtV -d\v{Vu) + [2ii\Du\'^ + \{dwuf], 



W := T^kdkW-V{u-Vu)-V 



(/iAu + (A + /i) Vdiv u) 

1 + a ' 



aV{n + V) 



(1 + a) 

.1 .7 

According to Proposition 14.11 we thus have to bound A, R in L^(i?2 i H -B^ i) ^'^'^ -^' ^ ^'^ 
L^{B2i n ^2 i)- We shall assume throughout that ||a||ioo™+ xr^-j is small. 

Bounds for \\A\\ i ? . Recall that 

^ = —Tq^^u — adivii. 
Using standard product laws for the paraproduct and remainder (see e.g. [2]), we get 
(4.93) 

(4.94) ||r^,„n«|| .1 <||Va|| .i \\u 



\TkaA 



1 < ||Va|| _i ||u|| 5 , 



(4.95) 
(4.96) 

Hence 

(4.97) 



\TkaU'\ 

lladivul 



< Hall 1 ||divM|| 3 , 



Li(45 ) ~ " "l-(b|i) 



1^1 



< 



L\B;\nBl^) 



I 7 + ||a|| 7 ||divn|| 3 . 

Li(B|,) i°°(B|i) iH4^i) 



1 7 m 5 9 . 



Bounds for \\D\\ 1 5 and ||Ty|| 1 5 . We may rewrite D as follows: 

L^iBl^nBl^) L^Bl.nBl,) 



D = [Tu„A-^d^]dku' - A-^diT'u^ - A-Miv 



1 + a 



(/iAu + (A + /i) Vdiv u-V{n + V)) 



The first two terms of D may be treated as in (j3.49p : we get for any s > 0, 
(4.98) \\[T^,,A-^di\dku' - A-^d,T'Q^^.u%^ < ||Vu|U^ ||n||^ 

Next, classical composition and tame estimates yield for s > 0, 



1 + a 



-^mI 



BS. 



< ||a||Loo||^n||o, + ||^M||Loo||a| 
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. 3 

Hence, using the embedding i?! i ^^ L°°, we easily get 

(4.99) II ^u|| 1 5 < ||a|| 1 5 Hull 5 9 . 

"1 + a \HBl,nBl,)-" "L^iBl.nBiy "LHBl^nBl,) 

Finally, we have 

(4.100) ||^v(7^ + y)|| ^ .1 <l|a|| .1 ||v(7^ + y)|| 3 

(4.101) ||,^v(7e + y)|| 5 <||a|| 3 ||v(7^ + y)|| 

+M .5 ||v(7^ + y)|| 3 . 

So putting ()i:98]) to (|4.10ip together, we get 

\\D\\ 1 5 < ||ii|| 5 ||tt|| 1 5 

LHBl.nBl,) L^Bl^) L--(Bl,nBl,) 

(4.102) , 

+ \\a\\ 1 5 (||ii|| 5 9 + 1|7^|| 5 7 + ||vy|| 3 5 

L--(Bl,nBl,y L^Bl.nBl,) LHBl.nBl,) L^Bl,nBl,y 

It is clear that W satisfies exactly the same inequality. 

Bounds for \\R\\ i 7 . Recall that 
L^Bi.nBi,) 

R = -T'q^t^u'' - 7^divn - dtV - d\N{Vu) + [2fi\Du\'^ + \{dwuf]. 

First we have for any s > 0, 

WT'd^n'^^WBi^ < l|V7^||Loo||u||^.^^. 

Hence 

iirr,; i-ii ^ 11,^ II II II 

\u\\ 



L^Bl,) L^Bl^) L--{Bl^) 



J-a^-RU II .1 >, II /-ell 



(4.103) 

II T^/ fell ^^llTl,, ,,,, 

11% T^n II .7 < ||7^|| .5 HUH .7 



Next, product estimates imply that 
(4.104) 



T^divull 1 < IIT^II 3 ||n|| 3 , 

T^divull 7 < ||7^|| 3 ||u|| 9 + ||u|| 5 ||7^|| 7 . 

LHBi,) L--{Bl^) LHBl,) LHBl,) L-o(bI^) 



We also have 



(4.105) ||div(yu)|| 1 < ||y|| 1 ||u|| 5 +||y|| 5 ||u|| 1 

L^Bl,) "Loo(Biy \HbI,) LHBiy \^{Bl,) 

(4.106) ||div(yu)|| 7 < ||F|| 3 ||u|| 9 +||y|| 9 ||m|| 3 . 

L^Bl,) L--(Bl^) LHBl,) LHBl,) L--{BI,) 

And finally, 

(4.107) ||Vu(g)Vn|| 1 < ||Vn||rwrcx)J|Vn|| 1 , 

(4.108) ||Vu(8)Vn|| 7 < ||Vn||roo/rooJ|Vn|| 7 . 
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Therefore, combining inequalities (J4.103P to (|4.108p . and using embedding, we end up with 



\R\\ 



< 



\u\ 



9 \\U\\ 



(4.109) 



5 . 9 I|m 



+ \\dtV\ 



LHBl.nBl,) 



+ \\V\\ 5 9 \\u\\ 1 3 + \\u\\ 5 9 \\V\\ 1 3 

LHBl.nBl,) L°-{Bl,nBl,) LHBl.nBl,) L--{Bl,nBl,) 

+ ||7^|| 1 7 ||li|| 5 9 +117^11 3 ||ti|| 3 +117^11 5 ||u|| 1 . 

L--{Bl,nBl,) LHBl^nBl,) L2(Bl,) L^bI,) L^bI,) L--{BI,) 



Putting ()4.97p . (|4.102p and (|4.109p together, one may finally conclude that for some con- 
stant K depending only on A and jl, we have 



X <k(x{o)+x^ + (\\v\ 



+ \\v\ 



)x + \\dtV\\ 1 



■L-(B|,nB2^i) L^Bl,nBl,y L^Bl^nBl, 

From this, we see that if X{0) and the terms pertaining to V are small enough, then 
(4.110) X <2K(x{0) + \\dtV 



LHBl.nBl,)^ 

Going back to the original variables according to (j4.88p , we then get the global existence part 
of Theorem 12.31 for any e > 0. 

4.3. The proof of convergence. As in the case where k > 0, we first show that {Qu^,TZ^) 
goes to 0, a consequence of Strichartz estimates, then establish that (Vu^,Q^) goes to the 
solution (f , 0) of the Boussinesq system (jl.lip . 

4.3.1. Convergence to for (Qu'^ ,TZ^). It suffices to prove dispersion estimates in the case 
e = 1. The change of variable (J4.88P will provide us with decay estimates in the general case. 
Now, the system for {Qu,TZ) reads 



dtQu + Vn = -Q{u- Vu) - Q 



2 /\ , ,,\/j-.„,\2 



dtJZ + div Qu = -dtV - div {{V + n)u) - 2fi\Du\^ - (A + fi){dwuy =: Ms. 
Therefore, Strichartz estimates imply that for all p £ [2,oo), 
(4.111) 



(Q7x,7^)|L^ .a_i <||(Qno,7^o)||.l + 



'■"■^'"..(4)- 



So it is only a matter of bounding Hi and EI2 in L^(i?|]^), which may be done by using standard 
results of continuity in Besov spaces and the fact that Q is an homogeneous multiplier of 
degree 0. More precisely, we have 



\Q 



Q{u ■ Vn 
Au 



l + a 



■ l + a 



\\dw{{v + n)u)\ 

||Vn(g) Vul 



< 



\u\ 



1 \\u\ 



LHB4^) 



L--{Bi,) L^{Bl,) 



LHBl,) 



< 



(1 + ||a|| 3 )||ti|| 5 , 

"l^iBl^" LHBi,) 



< 



LHBl,) ~ L-(i?2^i) 



\vv\ 



LHBl^) 



+ \\V'R\ 



LHBi.y 



< 



LHBl,)^" "l^bI,)' 



5 \\v + n\\ 



I + \\u\ 



L-°iBi,) L°-{Bi) 



1 \\V + Tl\\ 



LHBi.y 



LHBl,) 



< 



\u\ 



5 \\U\ 



LHBl,) L--{Bl,) 
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Therefore, if we set 

Co =||(ao,'^o)ll 1 7 +11^011 1 5 +115*^11 1 7 , 

then plugging the above inequahties in (j4.11ip and using (j4.110p leads to 
\\{Qu,n)\\ 2p 2_i <KCo for all pG [2,oo). 

. 5 

From (j4.110p . we also know that {Qu,lZ) is bounded by KCq in L^lB^i)- Hence using 
interpolation exactly as in the case k > leads to 

\\{Qu,lZ)\\i2, r.s ^<KCq for ah p> 2 and s G [-1/2 + 4/p, 3/p]. 
Now, going back to the original variables, we gather that for e > 0, we have 

(4.112) ||(Qu^7^^)|| 2p 2_i <KC0'^ if 2<|j<oo, 

L^{Bl^^) 

(4.113) l/5||(gu^7^")||f2fo. . < KCI^{ev)l~' for ah p>2 and s G [-1/2 + 4/p,3/p]. 

with Cq defined in the statement of Theorem 12.31 

4.3.2. Global existence of a solution to (jl.lip . Under the assumption that ()2.17p . the ex- 
istence of a global solution (0, v) to (jl.lip satisfying ()2.18p is an easy modification of the 
corresponding proof for the standard incompressible Navier-Stokes equations, combined with 
the following a priori estimate for the transport equation (see e.g. [2j, Chap. 3): 

l-T 

\\v\\ . 5 dt 

^T >.^2,1-' ^2,1 V-'O 

Indeed, using once again Proposition 15.21 and product estimates, we see that 



leii 1 < iieoll 1 expi / 

L-(B|i) Bl^ \Jo 



\\u\\ 1 +;U||n|| 5 ^ Ikoll i +11^11 1 IKII 5 + ||0|| 1 ||Vy|| 3 . 

Hence if (j2.17p is fulfilled then one may close the a priori estimates globally in time. 
4.3.3. Convergence of {Vu'^, 0^). Let us first notice that (recall that 0^ = a^ — TZ^ — V^) 

= p(0^v(y^ + 7^")) -pu^J^v(y^+7^" 

Therefore the system for {dd'^,dv'') := (0^ — Q,Vu^ — v) writes 
' at(50^ + Vu^ ■ Vdd^ = -dv^ ■ V0 - Qu^ ■ V0^ - 0^div Qu^ - e(2^|Z)n=|2 + A(divM=)2), 

dtdv^ - fiASv^ + V{v • V&") + V{dv^ • VVu") = V{m^VV + Q'VSV^ + 0^V7^^) 

-v(qu^ ■ VVu^ + u^ ■ VQu^ + j^{Au^ + a^ViV^+n^))]- 

In contrast with the heat-conducting case, we do not know how to prove convergence globally 
in time. This is due to the fact that some terms in the right-hand side of the equations 
for ((©^,(5u^) decay to only in Lp'-in time spaces and that (59^ satisfies a mere transport 
equation (hence the r.h.s. should be bounded in V- -in-time space if we want to get a time 
independent bound for (S9^). 
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We claim nevertheless that 0^ — )■ in L'^^{B't~^) with s as in the previous step, and that 
Vu^ — )■ 'u in 

Let us first examine ^^ . Denoting by Ki the r.h.s. of the equation for i©"^, standard estimates 
for the transport equation ensure that, if s > —1/2 then we have for all T > 0, 

(4.114) ||<©1lfoo.RS-2^ < expf / llVT^u^ll 3 dt ) ( ll*>nll r=-2 + / ||IKi|Us-2 dt ) • 

Product laws give if, in addition, s > 1/2, 

ll*'-ve||^.-2 < 11*^11^. j|VG||^_i, 

-°2,1 

||GMivQ'u^|U.-2 < lldivQ-u^lLs-ilie^ll i . 

For the last term of Ki , we use the fact that the product maps i?2 1 ^ ^2^ 1 ^^ ^2 1 ^^ 
< Q < 1. Hence using the embedding B21 M- -Bp^^ with a = 3/p — s, we get 

||2/.|Z?n^|2 + A(divn^)2||^... < ||n^|| , ||^-|| ,.. 

p,l R ^ R ^ 



Inserting those inequalities in (|4.114p and keeping in mind that VPu^ is uniformly bounded 

1+4 31 
2 ^ p' pJ 



in L^iB^^), we get for any s G [-i + f , f] n (^,00): 



|(5B^||foo^RS-2^ < ||«)oI|r«-2 + / ||*^||rs ||0||.i di 



2,1 

'0 V" "^P,i" "b|, "b^^/ "b|,"' 



whence 



(4.115) ii<©^ii^^(^._2. < \m\\^s-2 

T ^ p,l I p,l 

+ (l + Ti)||eo|| 1 ||&^||^,(^.^)+^w^.^) + z.-nC7o^)'((^T)i(ei.)° + eKz.T)t). 



In order to bound (fc^, we shall make use once again of the parabolic estimates given by 
Proposition 15.21 The main difficulty here is that some terms of the r.h.s. 1K2 of the equation 
for bv^ cannot be bounded in global L^-in-time spaces. Hence we shall use the following 
inequality which may be easily deduced from Proposition l5.2l (we do not track the dependency 
with respect to //): 
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Now, from product estimates in Besov spaces, we get 



v-Vdv^ 



^L'riKTi'+Kf) 



SU^ • VPU^ILl ^RS-l I RS-2N 

,B"Vy||RS-2 

Q'VSV'\\f2lf,s-2. 
QU" ■VVu'\\rl,f.s-l., 

"^t(^p,i ) 

u^ ■ VQu^Li /Rs-i-, 
^Ty-°p.i ) 



i\\v\ 



3 + \\V\ 
r2 (n^ 



< llVyll 3 ||(5B^i 



^L¥iB;:i'+K:i'y 






B^ 

-°2,1 



I - -1 



|V(5l/^| 



^?^(^ri')' 



L-{BJ,) 



ivTe" 



^L^Kri'y 



< llVPu^ 



^t(-^2,1 



IQ-u^ 



I^|(-B;,i)' 



< 



k'^ll _ 3 ||VQu^||r2 ^Rs-i-, 



■^t(-^2,i) 



li^c^ri )' 



and arguing as in the proof of ()3.72p , 



tAw 



l+ea--^'^ llLi,(B^-/) 



< 



ea' 



L^{Beu ) Lj,{B2i) 



W^u'W ^ .1 

-'^t(-^2,i) 
,£|| 



Finally, because -62,1^ ^^ -^p,i^' 



|^a-V(y^+7l-)||^^(^.-2)<||ea^lL ., 

<i/-i(ei/)°||a^ 






L^i^eJ ) L5^(B2i) 



1 F^ + 7^^ 



■^t('^2,i) 



Therefore, putting together all those inequalities and using the estimates provided by the 
previous steps we conclude that 



I*" 



i^T (s^i +^ri ) 



+ 11*' 



^lUb;,i)+lUbi,) - ii^^oiib;- 



< \\dVn\\t,s-l,t,s-2 + 



p,l 



T 







RS-2 



If v^^Cq is suitably small, we thus deduce that 

with /3(e) := ||&§ 11^.-1^^.-2 + Cg(||V,5V^^||^^(^.-i) + Cge"). 

Therefore, plugging ()4.115p in the above integral, and using Gronwall lemma, we get 

(4.116) 



I'»'IIl-(b-^) < ll«>Sb-2 + (i + T2)||eo|| 1 (IIHb-i+B-2 



1 3. 



+C7o^(||V<r^||^^(^.-i) + C^e-^-')) + {Cf,f{T-2ev-' + eT^^'^ ) 



xexpi^iieoll 1 (i + r2)||vy|| 3 



^2,1 
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and 



(4.117) 11*^11^^(^.-1^^.-2) + W'\\mB;^,)+L},{B;^,) 



p,l ' p,l T^ p,l ' 



+K ||«§||^.-2 + (1 + Ti)||eo|| 1 (iiHb-^+B-^ 



+c'g(l|V(5v^^||^2 (^.-1) + Q^e?^^)) + (Co^)^(rtep-^ + eri-t) ) ||vy|| 3 

P' / ^T\ 2,l) 

xexp(K||eo|| 1 (i + r5)||vy|| 3 

V 4^1 ^Hi^2^l) 

whenever s G [—1/2 + 4/p, 3/p] and s > 1/2. This ensures the convergence of (0^,'Pu'^) to 
(0,t;) with an exphcit rate. ■ 

5. Appendix 

In this Appendix, we give some a priori estimates involving hybrid Besov spaces. Let us 
start with product estimates. 

Lemma 5.1. Suppose that p £ [2, oo] and /3 > 0. There exists a constant C such that for all 
a > 0, we have 

Wfgh^-^^- < II/IIb-II5||.3 , if (3-1/2 <s<3/p, 

■^2,1 

11/511^-/3.+ < Il/lb-+ll5ll.3 , z/ /3-3/2<s<3/p-l. 

-°2,1 

Proof. We may assume that a = 1 making a change of variables if the case may be. In order 
to prove the first inequality, it suffices to notice that for all o" G M we have 

• ■ — — 3 • —R 

Now, it is well known (see e.g. [2]) that the usual product maps B"^ x B21 in B^^^ 
whenever /3 — 3/2 < a < 3/p and /3 > 0. Therefore 



II/5||r--/3-i < ll/llij--i||5'IL3-« and \\fg\\r,s-fi 



< 



is-/^-- r.. iuiIb--i||5'II .3 ,3 and 11/511^-/3 ;:; II Jlls« 11511 . 3 ^. 

p,i P.i B^-^ p,i P.i Bj^-^ 

This implies the first inequality. 

Proving the second inequality is rather similar: now we use the fact that 

(5.119) • Bs,+ ~ • RS _l_RS+l- 

Decomposing / into low and high frequencies according to ()2.12p . we have 
Now, the aforementioned product law ensures that 



3 

B^' 

Dry -, 



Wf^gWs^+i-P < ||/b=+i II5II . 3 ;3 and \\f^g\\^s-0 < Wf^W^s \\g 

p,i P.i Bj J P>i P'^ 

So taking advantage of ()5.119p completes the proof of the second inequality. ■ 

The following Strichartz estimates for the acoustic wave equation are the key to the proof 
of convergence. 
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Proposition 5.1. Let (q, Qu) (with cuilQu = 0) satisfy the 3D acoustic wave equation 

j dtq + V2divQu = F, 

\ dtQu + y/2Vq = G. 
Then for any a > 0, s G M and p G [2, oo) the following estimates hold true 

\\iq,Qu)Lj, ..+!_, <c(\\{qo,Quo)\\^s+\\iF,G)\\^,^^^X 

^ p,l ' 



Proof. The first inequality has been proved in [6|. In order to prove the second one, one just 
has to decompose {q, Qu) into low and high frequencies, that is {q, Qu) = (g^, Qu^) + {q^ , Qu^) 
and apply the first inequality with s it 1 (resp. s) to (g^, Qu^) (resp. {q^, Qu^)). ■ 

Let us finally state maximal regularity estimates for the heat equation, in hybrid Besov 
spaces. 

Proposition 5.2. Let u he a solution to the heat equation 

dtu - Au = f, 

u\t=0 = Uq. 

Then we have the following estimates for any a £ M, a > 0, p G [1, oo] and q > r: 



\\^\\~ .<r+^ ^ lko||^<^ +IUII~ .<T+^-2,l 

Proof. The first inequality is classical (see e.g. [2], Chap. 3). The second inequality may be 
obtained from the first one after decomposing u, uq and / into low and high frequencies. ■ 
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